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Abstract 


This report details the computer hnplementation and calculations 
of the induced velocities produced by a wake model consisting of a 
trailing vortex system defined from a prescribed time averaged down- 
wash distribution. Induced velocities are computed by approximating 
each spiral turn by a pair of large straight vortex segments positioned 
at critical points relative to where the induced velocity is required. 
A remainder term for the rest of the spiral turn is added. This ap- 
proach results in decreased computation time compared to classical 
models where each spiral turn is broken down in small straight vortex 
segments. The model includes featmes such as harmonic variation 
of circulation, down wash outside of the blade and/or outside the tip 
path plane, blade bound vorticity induced velocity with harmonic 
variation of circulation and time averaging. The influence of various 
options and parameters on the results are investigated and results are 
compared to experimental field measurements with which, a reason- 
able agreement is obtained. The capabilities of the model as well as 
its extension possibilities are studied. An assessment of the model for 
possible adptation to interaction problems shows it not to be suitable 
for this purpose. However, further features like a shed vortex sys- 
tem, an inboard trailing wake system, close blade vortex interaction 
and lifting surface effects can be added to it. The performance of the 
model in predicting the recently-acquired NASA Lcingley Inflow data 
base for a four-bladed rotor is compared to that of the Scully Free 
Wake code, a well-established program which requires much greater 
computational resources. It is foimd that the two codes predict the 
experimental data with essentially the same accuracy, and show the 
same trends. The model is an efficient and flexible tool to obtciin a de- 
scription of a rotor’s induced downwash field in a workstation- type of 
environment with accuracies comparable to those obtained by codes 
requiring an order of magnitude more of computer resources. 
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1 Nomenclature 


aoc- Zeroth harmonic coefficient for the bound circxilation 
aic". First harmonic, longitudinal, for the bound circulation 
bic’ First harmonic, lateral, for the bound circulation 
b: Number of blades 
c: Blade chord 

c': Blade chord, normalized by R 
Cc‘. Vortex core radius, normalized by c 


-ro/liric V\Tr ?? ( 

^ ^ A4.VP/4. J \ ^ JX 


r.'\ 
' / 


Ct’> Thrust coefficient T/(/07ri2^fl^) 


d: Weighting coefficient relative to remainder term contribution to 
induced downwash computation 


*Wa<fc' Blade index ( range: 1 to 5 ) 


Spiral turn index ( range: 1 to iV ) 

i: unit vector parallel to tip path plane and to the vertical plane 
containing the freestream velocity 

i^c •= (dOm/dAV’v)/||dOTn/dA^„|| unit tangent vector to the spiral 
at a given point, 

j: unit vector parallel to tip path plane and such that j = k x i. 

jAC- unit normal vector to the spiral at a given point, {jac = kx Iac) 

k: unit vector perpendicular to the tip path plane and directed up- 
ward. 


N: Number of spiral turns used for the trailing vortex system of each 
blade. 
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• ^MEAN‘ Number of divisions of one sector used to average the in- 
duced downwash over azimuth. (One sector is the angle between two 
consecutive blades, i.e. (27 t/6) 

• P; Point where the induced velocity is being computed 

• pb: projection of vector PB on jac (pb = P~BJac) 

• Tp: Radial coordinate of point P, where the induced velocity is being 
computed normalized by R 

• r„: Contracted radius, normalized by R 

• R: Rotor radius 

• R: Algebraic curvature radius of projected spiral in tip path plane 

• T: Rotor thrust 

• v: Induced velocity (positive up) 

• V: Free stream velocity 

• X, y, z: Coordinates of a point in the i, j, k vector base, 

• zp: z ordinate difference between point P and tip path plane, nor- 
malized by R, (Positive above tip path plane). 

• zpB‘ z ordinate difference between point P and point B, normalized 
by R {zpB = k.PB/R = zb - zp = z^ — zp) 

• a: Angle of attack, of rotor tip path plane (negative in forward flight) 

• a: Radial abscissa integration variable 

• /?: Angle between vortex segment and blade of reference {/3 = ^p) 

• f3: Integration variable 

• A^t,: age of a vortex hne point. 
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• An Radial length used to average the induced downvpash when such 
averaging is used normalized by R. 

• F: Circulation strength 

• Fj,: Circtdation strength of index ibiade at azimuth ipb 

• F': Normalized circidation strength (F/(i?^fl)) 

• Fe//: Corrected circulation strength with respect to vortex core model 

• k: Wake skew angle, (— tan"^(/ij,/A)) 

• A: Inflow ratio, + Aj), positive up. 

• Aj: Induced inflow ratio, (v/fIR) 

• A*: Reduced induced inflow ratio, (Aj/AjoH) 

• Ajo: Momentum value of induced inflow ratio 

• AjoH’ Momentum value of induced inflow ratio at hover ( — y(7r/2) 

• Hx- Advance ratio, (F cosa/fii?) in tip path plane reference frame 

• Complementary ratio, {V sin a /UR) in tip path plane reference 
frame, negative in forward flight or chmb. 

• fl: Rotor rotational velocity 

• Azimuth angle of straight vortex line segment ($„ = axgMc) 

• Azimuth angle of blade of index iuade 

• ij>p: Azimuth angle of point P where the induced velocity is computed 

• Azimuth angle of reference blade 

• Azimuth angle at which a vortex line point was shed by a blade. 

• p: Radial distance from vortex center at point B to point P, divided 
by core radius 

R 
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• p: Air density, where appropriate 

• <r: Rotor solidity, {hc/irR) 
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Figure 1: Induced dovmwash computation 

2 Introduction 

Beddoes [1] showed that the problem of determining the induced downwash 
through a rotor at low advance ratios can be solved in two successive steps. 
First, a geometric model of the distorted trailing vortex system of the wake 
is defined and computed from a prescribed time averaged downwash dis- 
tribution. Then, the induced velocity is evaluated by positioning simple 
straight vortex segments at’ critical locations in the wake with respect io 
the point where the induced velocity is computed and adding a remainder 
term to account for the neglected portions of the wake. The computation 
time of the induced velocity at a given point can thus be divided in two 
parts: Location of the critical positions of the wake (two for each spiral 
turn) and computation of the three terms making up the induced velocity 
for each spiral turn (a straight vortex segment for each of the two posi- 
tions and a remainder term). The advantage of this approach is decreased 
computation time compared to a classical scheme where each spiral turn is 
broken down in small straight vortex segments, the location of which must 
also be computed. 

An induced velocity computation is characterized by the point P where 
the induced downwash is required and the rotor position defined by the 
instantaneous azimuth of a reference blade V’r (Figure 1). P can be placed 
on a blade to compute the instantaneous downwash seen by the blade at 
a given azimuth, or outside the blade and even outside the tip path plane 
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Figtire 2: Trailing vortex system 


3 General description of the wake model 

The wake modeled here is the trailing vortex system only. A reasonable 
agreement with experimental results is obtained without a shed wake which 
would be resource-and time-consuming. The trailing wake is represented 
by effectively rolled up tip vortex lines trailed by each blade from effective 
vortex radial sources at location r„, contracted radius. Each of these spiral 
Knes foHows a cyUnder skewed according to the advance ratio (Figure 2). 
The model defines the geometry of those lines and then computes the in- 
duced velocities by approximating each spiral turn by a pair of large straight 
vortex segments positioned at critical points relative to where the induced 
velocity is required. A remainder term for the rest of the spiral turn is then 
axlded. It is derived from a closed form expression of the velocity induced 
by a circular vortex ring. A finite core model is used for the computations. 
Radial averaging of the computed downwash over a sampling radial length 
along a blade can be used to remove peaks due to arbitrary sampling. Those 
peaks occur when a vortex segment happens to be positioned very near the 
point where the induced velocity is being computed. It is then considered 
to be a numerically induced peak which does not necessarily have a signif- 
icant physical meaning because it is an occurence related to the particuleu: 
model employed and not to the passage of a vortex line at that particular 
point. 


UiXlEDlNG PAGE BLANK NOT FILMED 
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4 Tip vortex geometry 

A prescribed geometry presented in reference [2] is used. The vertical dis- 
placement of elements of the wake is computed by integration of a pre- 
scribed time averaged downwash field given in reference [3]. The tip vortex 
lines are represented parametrically in the tip path plane frame of reference 
(i, j, k) moving with the rotor. The parameters are: 

1. hiade or the index of the blade from which the given vortex line trails 
( range: 1 to 6 ) 

2. the angular position of the reference rotor blade (iuotU = 1) ( range: 
0 to 27 t ) 

3. AV’t, the age of a vortex fine element. It is defined as the angle by 
which the rotor has turned since that element was produced. ( range: 
0 to 2N7T ) 


Thus, 

Xv = r„ cos + /i* AV>v 
^ yv = rv sin V’t, > 

. A# J 

with 

V’w — V’r "1“ ‘^'^{^blade 1)/^ 

It follows from the definition X = fiz + that 

+ / Xidij) 
Jo 


( 1 ) 

( 2 ) 

( 3 ) 


The averaged downwash field, within the disc defined by the contracted 
radius is given by: 


Ai = A.o(l + Ex^ - E\y'^\) 


( 4 ) 


where primed values are normalized by and E = l«l* Outside the disc, 
and dowstream of the tip path plane center, the averaged downwash field 
is given by: 


Ai = 2Aio(l - E\y% 


( 5 ) 
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Here, A^o is the momentum theory value of the induced inflow ratio and is 
the solution of: 

Three cases present themselves: 

1. cos^„ > 0: the vortex element has always been outside the disc, and 
downstream of the tip path plane center (Figure 3). 

2 „ = fizAip„ + 2Ai, / (1 - E\y'^\)dip 

Jo 

2 „ = + 2Ai„(l - E\y'^\)A-tf)„ (7) 

2. cos < 0, «„ < — r„ cos the vortex element has not moved beyond 

the disc (Figure 4). 

2 „ = fizA-ip^ + ko (1 + Ex' - E\y'^\)dt{} 

Jo 

+ Xi^ / (1 + E{cosrJ)^ + - E\y^^\)d^ 

Jo 

2 „ = /x^AV>t, + Ai„(l + E{cosip^ + /ij.AV’v/(2r„) - |t/'®|))A^„ (8) 
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3. cos^„ < 0, x„ > —r^cos'ip^'. the vortex element has spent part of 
the time within the disc and the remainder outside of the disc, and 
dowstream of the tip path plane center (Figure 5). 

/ {l + Ex'-E\y'^\)di> + 2XiJ (1 - E\y'^\)diP 

Jo JA\I>vi 

where AV>vi = cos^„//ij, 

Thus: 


+Aio(l +^(cosV>„ +/i*(-27'„cosV’v/A‘x)/(2r„) - |y'®|) 

(- 2 t’„ cosV>„/Mx) + 2 Ai„(l - f;|y'®l)(AV’« + 2r„ cos V’„//iiz) 

which simplifies to: 


Zy — 2AtQ(l E^\y |)®v//^x 


( 9 ) 
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5 Computation of induced velocities 

For each spiral turn, the vortex hnes trailed by each blade are accounted 
for in the computation of induced velocity by placing a pair of straight 
vortex segments at critical locations, and adding a remainder term. In this 
model, the vortex lines have first harmonic circulation variation with the 
local a.zimuth V’t> at which a particular point was shed. A point belonging 
to a vortex line will be attributed a circulation strength equal to that of 
the blade from which the line is trailing at the time when the blade shed 
that paticular point, that is azimuth The zeroth harmonic component 
of circulation is related to the thrust coefficient. Since each blade is trail- 
ing concentrated vortex lines from a single radial station (the contracted 
radius), then the bound vorticity should be constant along each blade sind 
then, fall off abruptly at the contracted radius. Since no inboard trailing 
vortex fine is included, the bound vorticity should be constant from the 
center of the rotor to the contracted radius. Thus, the bound vorticity of 
each blade is modeled as constant circulation straight vortex segments of 
length equal to the contracted radius. 

Reference [1] presents a radially averaged downwash expression. It can 
be used for the distribution of downwash along a blade to smooth out peaks 
appearing because of the discrete nature of the sample of points where the 
induced downwash is computed. Those peaks occur when a vortex segment 
happens to be positioned very near the point where the induced velocity 
is being computed. It is then considered to be a numerically induced peak 
which does not necessarily have a significant physical meaning because it 
is an occurence related to the particular model employed and not to the 
passage of a vortex fine at that particular point. The averaging length then 
shotild be of the order of the definition interval of the sampling. 

5.1 Pair of vortex segments 

For a given point P where the induced velocity is needed, and for each 
spiral turn of each vortex fine, two points B axe positioned to support a 
straight vortex segment of length 2r„ and circulation T^ff. The segments 
are horizontal (i.e. oriented parallel to the tip path plane) and each is 
contained in a vertical plane tangent to the spiral such that the vertical 


5 COMPUTATION OF INDUCED VELOCITIES 


17 



orthogonal plane passing through the point of tangency contains the point 
P (Figure 6). Thus, the trace of the second plane in the tip path plane 
joins the projection of point P to the point of tangency of the projected 
spiral. The velocity induced by say, segment AC, placed at B, at point P 
can be computed as follows: 

fM=c 

/ =: ( 10 ) 


r ^eff 
VB = k — — 
47T 


lIMPf 

where k is the unit vector perpendicular to the tip path plane and directed 
upward, dl = dlij^c — Rdaijic^ where is the unit tangent vector to the 
projected spiral at the projection of B in the tip path plane: 


dOm 

where Om — x^i + t/„j. From equation (1), we obtain: 

dOm , . , sv . , 

-jT-r- = (r„ sin ipv + fix P + { cos ip^)j 

Defining j^c ^ normal vector to the projected spiral at the projec- 

tion of B on the tip path plane: 

dlAC 


( 11 ) 
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Figure 7: Local unit vectors 

Jac is by construction directed toward the outside of the concavity of the 
projected spiral, at the projection of B in the tip path plane (Figure 7). For 
advance ratios smaller than 1, there is no inflexion point on the projected 
spiral and thus jjic is also: 


JAC = k X ij^c 


Fe//., in algebraic value around is gi-ven in the form based on refer- 
ence [4], 


r.„W = r(:^) 


and 


r = ro(l — aicCosV>v — feicsin^v) 
where p is the ratio between radial distance and core radius: 

2 (Pb + ^pb) 

P = 


( 12 ) 

(13) 

(14) 


Reference [1] suggests a constant circulation strength given here in non- 
dimensional form: 

V - ^ . . 
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Figure 8: Vector decomposition 


with Oo = 2.4. Indeed, the total thrust of a rotor having constant circrila- 
tion, full length blades could be expressed as follows: 


T = bpT I Drdr = bpTQ, 
Jo 


B? 

2 


Thus, 


which yields: 


Ct = 


bT 


pirR^iDRY • 2TT.R^n 


r = 


R^H b 


(16) 


Looking at equation (15), we see that the value recommended a© = 2.4 ac- 
counts for losses of blade root cut out, tip loss and non constant circulation 
compared to the ideal case of equation 16, which corresponds to a value 
0 © = 2 . 

Now, projecting MP on the Jac^ k vector base, (Figure 8), 

MP = MB — Rzpgk — RpbJac 

where zpB = k.PB/R = zb — zp = z^, — zp, pB = PB.jac 
Then, 

MP = —RaiAC ~ RzpBk — RpbJac 
dlxMP = dlR^zpB] AC ~ Psk) 
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and 

k{dlxMP) = -Rpsdl = -R^peda 
Also, we have (MPy = (a^ + Zpg + p%)R^ thus: 

Te/y /■“=+’•«' pBda 

4xi?yo=-r. + 

Switching to non-dimensional variables with Vg = vb/{RD) and T'^ff = 
Teff/{R^Cl)^ we obtain: 


Vt) — 




ya=+r,, 

Ja=-rv (CK 


47T 


da 


(a2 4 b +Pb)®^^ 


Letting yS* = 4 b + Pb sinh^ = 


, _ KffPB /•“■Bsinh(^) d9 

Jo 


Vb = -- 


2ir^^ 

, _ _£^jtP£ 


cosh^ 6 


Vg = 


27 t /?2 tanh(argsinh(^)) 


t?D 


^eZ/PB 


1 


27T/?2 + 

Thus, we finally have the expression v'g: 

• ^'efj‘’'vPB 1 


Vg — 


(17) 


2'fr(zpB +Pb) (^bb + Pb + 

with r;„ = r(^) 

5.2 Remainder term 

The remainder term accounts for the contribution of the remaining por- 
tion of the spiral turn. Reference 1 provides the following approximate 
expression. 


”x r,y,r, ^ ^ (1 _ pj/r,)jp^)(l + zJ.p)'/2 


(18) 
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with 

Now, 


Fft = ro(l — Oic cos - hic sin 


( 20 ) 


MP = MO + OP = —Raioc + R{xpi + ypj + zpk) 
and d/ = Rdaioc with ioc = cos + sinV’bJ. Thus, 

k{dlxMP) = R^{yp cos — xpsmrl>b)da 


and 

So, 


{MPf = ((xp - acosV>i,)^ + (yp - asinV>(,)^ + ^p)R^ 


V 

‘Wait 


^a=r» 

Ja=0 


(yp COS — xp sin V>i>) x 
da 


((xp — a cos^b)^ + (yp — a sin + Zp)V^ 


( 21 ) 


Now, let 


X(a) = (xp - acosV>b)^ 4- (yp - asin^i,)^ + z|, 

Developing, we have 

X (q) = — 2(xp cos V»b + yp sin + (®p + yj> + Zp) = A + Ba + Ca^ 

for which notation, an expression of the above integral can be found as: 

-4t) 


ra=r, do, 2 . 2Cr„ + B B 

,=o “ e' " vC4‘ 


with Q — 4^ — Thus, finally, 

, r;,/., , . ,,2,2Cr. + B B, 

= — (s,PCOs^s-xps.n,fs)g(-^ - 


( 22 ) 


with 


Q=AAC-B^ 


A = Xp + yp-\- Zp 
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B = -2(apcosV»6 + yp sinV>b) 

C = 1 

X[r,)=A + Br,+erl 

also, = r'(^) and, 


^2 ^ (OP X ioc)' 




2 _ [(gpt + ypj + zpk) X {cosipbi + sinV>b;)]^ 
\{xp sin V’i, - t/p cos V>b)fc + 2p cos - zp sin 


2 [(®p sin ij}b - yp cos V’b)^ + 2 p 

^ = -r 


(23) 


5.4 Radial averaging of downwash 

The computation of the downwash along the blade involves occurences of 
close blade-vortex interactions. This may residt in spurious variations of 
downwash dependent on radial grid definition and offset. To address this 
situation, a radially averaged expression of the downwash induced by the 
pairs of straight vortex fine segments over a radial interval may be used. 

The downwash induced by a straight vortex line segment located at 
point B on point P situated on the blade is given from equation (17) as: 

, _ p> 1 

1 + 27 t{zIb + p%) (zj>s + p% + 


withP^/^ = r'(Y^^)andp 
Then, 


(pj+4») 


Let T> = •2pp+Cji and £ = Zp^+T^. 


rv„ 


Vp — 


PB 


2ir (pI + C)(p| + £)■/! 

For an angle (3 between the vortex segment and the blade, the downwash 
expression at the control point P is averaged over the projected length 
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Ar' = I Ar. sin /3\ where Ar is the averaging distance, and /3 = where 

is the azimuth angle of the straight vortex line segment, ($„ = argt^c) 
Thus: 

_ rV„ rPB+ArV2 pgdr 

27rAr' JpB-Ar'/2 {p% + V){p% + 5)^/2 
Reference [1] provides the expression resulting from this integral: 


V 


f 

B — 


rv. \ (pl. + ey^^-k jpl+ey^^ + k - 

2n2kAr' ” [{pI^ + 5)1/2 + + 5)1/2 _ k 


(26) 


where pp^ = Pb ~ Ar'/2, ps^ =Ps + Ar'/2, k = {S - Vy/^ = (r2 - c]yy/^, 
and V — zpg + c^, 5 = Zpg + r2. 
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6 Positioning of the vortex elements 

For each spiral turn trailed by each blade, the position of the two points B 
have to be known by the age of the vortex element they axe at so that ps, 
zpB and which are necessary for the computation of induced velocities 
can be computed. 

The original method involves an empirical formula used for the location 
of the critical elements of the trailing vortex system. However, its im- 
plementation involves the use of an extra remainder term for the induced 
velocity computation as specified in reference [1] due to a discontinuity in 
the formulation. Its implementation also revealed it to be unsuitable for 
points where the induced velocity is wanted outside the blade and/or out- 
side the tip path plane. The method chosen for its simplicity and absolute 
reliability in all cases is an iteration method by dichotomy of each spiral 
turn where the two points B are to be found. The criterion for convergence 
is expressed in terms of the angle between the projection on the tip path 
plane of PB and Jac which, in absolute value is prescribed to be smaiUer 
than 3 degrees. It is judged that this value is sufficiently small to assure 
accuracy and repeatability of the computations, yet not so small that it 
would result in too many iterations and thus computation time. The al- 
gorithm first advances along the spiral turn in increments of sixty degrees 
until the above angle changes sign. Then, on the interval thus located, 
a dichotomy process is performed until the tolerance in angle is achieved. 
The other point B is then found on the remaining portion of the spiral 
turn in the same way. Let a step be defined as a computation of the angle 
between the projections of PB and Jaci that is, either after an increment 
along the spiral turn or in a division of the interval where there is a change 
of sign. With the increment of sixty degrees chosen, the number of steps 
necessary for convergence for each point has been minimized to six, in av- 
erage. The advantage of this method is a simplified algorithm, with an 
adjustable tolerance of positionning of the critical elements, at the same 
amount of computation time as the original method. No discontinuity ef- 
fects axe brought in by the algorithm, and thus the need for the extra 
remainder term is ehminated. However, the main reason for developing a 
new algorithm is that the original one could not handle points outside the 
blade, even after extensive modification. 
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7 Development of the code 

7.1 Structure of the code 


Care has been taken to separate the computation core of the code from the 
pre-and post processing of the computation. In addition, the modidarity 
of the wake model methodology was respected by separating code mod- 
ules dealing with the computation of the critical locations in the trailing 
vortex system from those dealing with the induced velocity computation. 
Thus, three core modules were created, each compilable, testable indepen- 
dently and from which, validating graphics can be generated: The first 
module finds the critical locations in the tip path plane of the trailing vor- 
tex system. The second module contains the subroutine for computing the 
momentum theory values of induced downwash as weU as the vertical dis- 
torted geometry of the trailing vortex system. The third module computes 
the induced velocity at a given point. The first two modules were tested 
through one program each, generating graphics shown in Figure 10 and 11, 
respectively. The third module which involves computation of induced ve- 
locity components is numerically tested through a third program with the 
debugger of the compiler. The computations axe traced out by extensive 
dumping throughout execution of variable and parameter values. These 
values axe then matched with hand calculations and correlation with geo- 
metrical insight available through the graphics generated by the first two 
modules. In addition, this last module is used by one front-end program 
producing the numerical results and generating the plots associated with 
those results. This program computes the induced downwash at one point 
of the flow field with or without averaging through the revolution of the 
rotor: A value Nmean of 1 corresponds to an induced downwash computa- 
tion at one point of the flow field with the axis of alignment of the reference 
blade passing through that blade. It is thus the instantaneous value of the 
induced downwash seen on the reference blade. A value Nmean of 2 cor- 
responds to the previous induced downwash computation averaged with 
another one done after the rotor has revolved one half sector. Experimen- 
tal data can thus be reproduced by choosing an appropriate Nmean (5 for 
example). This program produces both three dimensional (Figure 12, 13 
and 14) and contour plots (Figure 15 etc. ) of the results. All programs 
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use the same input file and can be linked and run together or at various 
degrees of independence. 

7.2 FORTRAN implementation 

The language used for coding is FORTRAN 77 with the exception of vari- 
able names. These use the capability of VAX/FORTRAN to recognize 
names of more than seven characters long. It is believed that this feature is 
of interest. Indeed, the variable names used in the three core modules each 
correspond to variables fisted in the nomenclature of this report. Thus, 
the use of names in the code of more than seven characters long makes it 
possible to make that correspondence immmediately, providing improved 
ease of reading the code. 

7.3 Graphics implementation 

As stated previously, the computation modules are completely independant 
for the pre and post processing. The post processing contains graphics- 
generating routines which make use of an elaborate integrated graphics 
package called DISSPLA ([6]). 
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8 Computations 

A test case for which experimental results axe available, (reference [5]), 
was treated. It involves a four bladed rotor, at an angle of attack of a = 
—3 degrees for an advance ratio ju* = .15. From this, a complementary 
ratio Hz = V sin a /CiR was computed. The input data indeed requires the 
specification of both Hx — ^ cos a/fliZ and fiz = V sina/QR a^ opposed to 
H = V/^R. 


8.1 Input Data 

The following input data are used for the three modules. The rotor charac- 
teristics, flight conditions and measurement plane location zp were chosen 
to replicate the test case of reference [5]: 

• Rotor characteristics: 

0 6 = 4 

o = 1. Contracted radius, normalized by R 
o cp = .07 Vortex core radius, normalized by R 

• Flight conditions: 

o Ct = .0064 
o Hx = .15 
o Hz = —.0078 

• Plane where induced downwash is computed: 

o Zp = .077 

Additional input data are necessary to cahbrate the model: 

• Xt, Contracted radius, normalized by R 

• cr Vortex core radius, normalized by R 

• N Number of spiral turns used for the trailing vortex system of each 
blade. 
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• Ooc Zeroth harmonic coefficient for the bound circulation 

• aic First harmonic, longitudinal, for the bound circulation 

• bic First harmonic, lateral, for the bound circulation 

• d Weighting coefficient relative to remainder term contribution to 
induced down wash computation 

• I AVERAGE Option Control integer relative to wether or not radial av- 
eraging computation of induced downwash is used 

• Nmean Number of divisions of one sector used to average the in- 
duced downwash over azimuth. (One sector is the angle between two 
successive blades, i.e. ( 27 t/ 6 ) 

• Ar Radial length used to average the induced downwash when such 
averaging is used 

8.2 First module visualization 

Figure 10 shows the critical locations in the trailing vortex system as seen 
in projection in the tip path plane for a given point P at •0P = 110 de- 
grees, Tp = 1.3 and zp = 0.077 and for a given rotor position = 200 
degrees. Each segment is indexed by iuodei and another index j = 1,2, 
respectively. There are two segments (j = 1,2) per spiral turn. 

8.3 Second module visualization 

Figure 11 shows the geometry of the distorted trailing vortex system as 
seen in perspective from downstream. 

8.4 Induced downwash computations 


Induced velocities are given in terms of A*, reduced induced inflow ratio, 

A* = Xi/XioH 
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where Aj is the induced inflow ratio: 

Xi = v/ClR 

aJid XioH is the momentum value of induced inflow ratio at hover: 

XioH = —\/Ct/2 

Thus, a positive value correspon ds to a n induced downwash and a value of 
1 corresponds to A,- = Xion = —\JCxj2 = —0.056 

8.4.1 Three dimensional plots 

Plots 12, 13 and 14 are three dimensional plots where positive values at 
a surface point mean a downwash. The parameters used axe the baseline 
values and they are modified in the subsequent plots. These parameters 
appear on the right-hand corner of each plot. 

• Figure 12 shows the plot of the induced downwash seen by a rotating 
blade from eizimuth V’r of 0 to 360 degrees along the span from .2 to 
1.2 radial position (non dimensionalized by R). The blade is seen to 
encounter regions of downwash everywhere except at the leading edge 
of the rotor disc in addition to ridges encountered frontally. 

• Figure 13 shows a plot of the induced downwash field on the region of 
azimuth 0 to 360 degrees along a radial position from .2 to 1.2 (non 
dimensionalized by R). Five averaging divisions of a rotor sector 
(Nmean) have been used. It can be seen that this plot is very similar 
to the preceding plot but it is smoother. Thus, the instantaneous 
induced downwash seen by a rotating blade and induced downwash 
field are very similar. This property is used in the subsequent compu- 
tations to evaluate the effect of the various parsimeters without going 
through the averaging process. 

• Figure 14 shows a plot of the induced downwash seen by a rotating 
blade from azimuth V’r of 90 to 180 degrees along the span from ,2R 
to 1.2il. This more detailed grid shows ridges encountered by the 
blade in a perpendicular fashion. 
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8.4.2 Contour plots, influence of input parameters 

The following plots are contour plots where positive values at a point on 
the rotor disc mean a downwash. Most computation meshes extend from 
,2R to 1.2R and the plots are truncated at I.IR. The lower and higher 
bounds of the radial positions of computed points can be changed at will, 
as can the truncation radial position for plotting. 

• Figure 15 shows the contour plot of the induced downwash seen by a 
rotating blade from azimuth of 0 to 360 degrees along the span from 
.2R to 1.1 J2. This plot contains the same information as Figure 12. No 
radial averaging is used, no bound circulation contribution is present 
and the first harmonic coefficients of the circulation are set to zero. 

• Figure 16 shows the contour plot of the induced downwash field on 
the region of azimuth 0 to 360 degrees from .2R to 1.1 jR. This plot 
contains the same information as Figure 13 and has the same param- 
eter settings as in the preceding plot except for Nmean which is set 
equal to 5, to average the velocities seen around the rotor as the rotor 
rotates by one sector in 5 increments. It is apparent that a region 
of upflow appears at the leading edge of the rotor disc. Compared 
with the experimental results obtained in reference [5] illustrated in 
Figure 34, there is a good agreement in the trends from fore to aft of 
the rotor disc. Also, an upflow on each side of the disc, downstream 
of 90 degrees and 270 degrees is obtained. The magnitudes of the 
downflow on the aft portion of the rotor are the same, providing an 
a priori confirmation for the values of oqc and d as valid. 

• Figure 17 shows that the application of the option of radially averaged 
induced velocity gives similar results than when not used (comparison 
with Figure 15). 

• Figure 18 and 19 show the influence of the number of spiral turns 
taken into account for the trailing vortex system. JV = 1 in Figure 18 
is too low. Neither trends nor magnitudes are approximating exper- 
imental results of Figure 34. However, it is seen that N = 5 gives 
almost the same results as JV = 3 (Figure 15). It is safe as a conse- 
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quence to retain three turns trailed by each blade for the subsequent 
computations. 

• Figure 20 shows that the contracted radius r„ is not a sensitive value 
for the induced velocity in the rotor disc by comparison with Fig- 
ure 15. 

• Figure 21 shows slightly more definition of the induced downwash 
seen by a rotating blade for a smaller vortex core radius cr although 
the sensitivity is not very high (comparison with Figure 15). It is, 
however the only parameter of the present model through which the 
aspect ratio effects of the rotor blades can be represented since the 
vortex core radius certainly bears some relationship to the relative 
size of the chord c and the blade radius R. 

• Figure 22 shows the influence of the flight condition variable fiz on the 
induced downwash seen by a rotating blade by comparison with the 
baseline case Figure 15. Zero freestream inflow component {/j,^ = 0. 
i.e ot = 0) seems, as expected, to bring more irregularities in the 
induced downwash contour patterns, since the trailing vortex system 
stays in the tip path plane longer than for a negative value of fix. 

• Figure 23 and 24 show the influence of the remainder part of the spiral 
turns contributions through the parameter d. It is seen to affect in 
a uniform way the whole rotor disc and the magnitude of .275 is 
confirmed as adequate. 

• Figure 25 and 26 show the influence of the zeroth harmonic compo- 
nent of the circulation through the parameter ooc- Here, it is seen 
that it influences the steepness of the fore to aft slope of induced 
downwash. The maximum downwash magnitude can thus be taken 
to size this parameter and the value of 2.4 is adequate. 

• Figure 27 and 28 show the influence of the first harmonic components 
of the circulation through the parameters au and b^. The asymmetry 
of the induced downwash pattern is clearly apparent here. 
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• Figure 29 shows the appearance of a larger region of upwash extend- 
ing from 160 to 270 degrees due to the value of aic in this induced 
downwash field plot. 

• Figure 30 illustrates the combined influence of Oic and 6i<, in this 
induced downwash field plot. The asymmetry, however does not ap- 
proximate experimental resiilts in a satisfactory way. 

• Figure 31 represents a possible choice of values of and bu giving 
a satisfactory approximation of trends and magnitude of the experi- 
mental residts obtained from reference [5] illustrated in Figure 34. 

• Figure 32 shows the small influence of blade bound circulation on the 
computation even with first harmonic variation of the circulation. 

8.4.3 Final results 

First harmonic coefficients of bound circulation were extracted from numer- 
ical data obtained with a free wake analysis code ([4]). Figure 33 shows the 
results obtained with those coefficients Oic = 0.25 and bu = —0.14. It is to 
be compared with experimental results obtained in reference [5] illustrated 
in Figure 34. An agreement in both patterns and magnitudes is clearly 
apparent: 

• A region of net upflow field appears at the leading edge crescent of 
the rotor disc. 

• A strong downflow exists only in the aft portion of the rotor disc. 

• An upflow on each side of the rotor disc, downstream of azimuth 90 
degrees and 270 degrees is obtained, 

• A certain degree of asymmetry is apparent with extremes in down 
and up flow on the advancing side of the disc. 

• The contour plots outboard of the rotor disc clearly show strong re- 
covery happening between 60 and 90 degrees on the advancing side 
and between 270 and 300 degrees on the retreating side which corre- 
lates well with experimental results. 
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• Finally, the magnitudes of the downflow in the aft portion of the rotor 
disc are the same as with experimental results. 

Figure 35 to 46 show comparison of those same resiilts with experimental 
data on azimuth lines every 30 degrees. The presence of a body and hub 
in the experimental setup in reference [5] might explain the larger region 
of upflow obtained in experimental results than for the computations. 

8.5 Computation time 

The computation of the induced velocity at a single point by the trailing 
vortex system takes one half second on a ’Digital VaxStation IF. The com- 
putation of the averaged induced velocity at a single point thus takes fives 
times this length (2.5 seconds). A field of 360 points will then be obtained 
in 15 minutes. 


8.6 Comparison with Scully Free- Wake Code 

For validation purposes, the Scully Free- Wake Code [4] was run for the 
NASA Langley test case. The radial variations of inflow velocity at different 
a.zimuths are shown in Figures 47 through 58. No significant trends are 
noticeable that differentiate the two codes. The agreement -with experiment 
is better at some azimuths with the Scully computation, and with the 
modified Beddoes method at others. At zero azimuth, the disagreement 
between the Scully code and experiment is physically reasonable: it may 
be argued that the disturbed flowfield downstream of the hub mast cannot 
be predicted by any code accurately. Surprisingly, no severe disagreement 
is visible for this case with the Beddoes method. At 300 and 330 degrees, 
on the other hand, the Beddoes method shows large error inboard, but the 
Scully code does not. At 90, 240, and 270 degrees, both codes perform very 
well. At most other azimuths, trends in difference with experiment are 
similar for both codes, yet the disagreement is not very large. It may thus 
be stated based on these comparisons that the accuracy of both methods 
in predicting the experimental test case is about the same. 
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9 Capabilities of the code 

9.1 Capabilities 

The test case computed with the present scheme (Figure 33) shows good 
agreement with experimental results obtained in reference [5] for a short 
computation time. Interesting results can otherwise be observed as to the 
similarity of instantaneous induced downwash seen by a rotating blade and 
induced downwash field plots. Different options and parameter values can 
also be used to study their influence on the results obtained and assess their 
sensitivity for a given rotor and flight conditions: 

• Utilization of radial averaging and the variation of the averaging in- 
terval used. 

• Variation of vortex core radius cr 

• Introduction of first harmonic variation of circulation strength over 
the trailing vortex system and on the blade bound vorticity. 

• Variation of the contracted radius r„ 

• Variation of the number of trailing vortex system turns taken into 
consideration 

• Variation of the number of increments over which the induced down- 
wash is averaged over azimuth (A reasonable number is 7). 

Factors to be considered in judging the completeness of the model are the 
presence of a shed vortex system, close blade vortex interaction, Ufting 
surface effects, inboard trailing wake system and presence of a body and 
hub as is the case for the experimental case in reference [5]. This last point 
might explain the larger region of upflow obtained in experimental results 
than for this model. 


9.2 Potential for growth 

The model implemented in this report is a simple wake model representing 
the trailing vortex system in an original way. Adding a shed vortex system. 
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an inboard trailing wake system, close blade vortex interaction and lifting 
surface effects, may not be consistent with the simplicity and peculiarity of 
the trailing vortex system model. The code might see its computation time 
increase to a level comparable to that of free wake codes and other much 
more refined wake models without having solid analytical consistency. On 
the other hand, adding those features might provide some insight on their 
relative importance in a wake model. 

9.3 Adaptability to interaction problems 

The wake model discussed in this report has its most fundamental charac- 
teristic in the geometry of its prescribed trailing vortex system. Each spiral 
trailing from the rotor has a projection in the tip path plane in the form 
of a regular cycloidal curve. This allows it to have the property of having 
only two tangents per turn which are perpendicular to the line joining the 
point of tangency to a given point where the induced velocity is computed. 
This property is fundamental to the present model consisting in computing 
the induced velocity by two staight vortex segments and a remainder term. 
It is also fundamental to the algorithm that searches for critical locations 
of those straight vortex segments. The interaction problem, consisting in 
distorting the shape of the trailing vortex system would result immediately 
in the loss of the geometrical property and render impossible the implemen- 
tation of both the model and its algorithm. Thus, this model appears to be 
unattractive for extension to problems involving complex vortex distortion 
effects. 
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10 Conclusions 

The test case computed with the present scheme (Figure 33) shows good 
agreement with experimental results obtained in reference [5] for a short 
computation time. The capabilities of the model include the study of cer- 
tain parameter influences on the results obtained and an assessment of 
their sensitivity for a given rotor and flight conditions. The model is by no 
means a complete wake model and it is possible to refine it and add some 
features Uke a shed vortex system, an inboard trailing wake system, close 
blade vortex interaction and lifting surface effects. Due to its peculiarity, 
the present model is not suitable to be evolved into a part of an interaction 
problem wake model. This model makes it possible to obtain a description 
of a rotor’s induced downwash field with very small computer resources, 
both in memory and time, in a workstation-type of environment. The per- 
formance of the model has been compared to that of the Scully Free-Wake 
code in predicting the extensive inflow velocity data base recently acquired 
at NASA Langley for a four-bladed rotor. The accuracy of the model is seen 
to be quite similar to that of the Scully code, and trends axe also similar. 
Yet the present model can be run on a workstation in interactive mode, 
whereas the Scully code requires a larger mainframe computer where it has 
to be run as a batch job for economic reasons. In conclusion, the model is 
an efficient and flexible tool for the preliminary evaluation of the induced 
velocity field of a rotor. 
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A User’s Manual 

A.l Introduction 

Three modules perform the computations. They are independent from the 
pre and post processing. The first module is Location and finds the critical 
locations in the tip path plane of the trailing vortex system. The second 
module is Wake and comprises the subroutine computing the momentum 
values of induced downwash as well as the vertical distorted geometry of 
the trailing vortex system. The third module is Velocity and computes 
the induced velocity at a given point. Each of the modules can be used 
through a pre and post processor: 

• Location is used by Test-Location 

• Wake is used by Test-Wake 

• Velocity is used by Test-Velocity 

Each module is thus testable independantly if the tests are effected in the 
order going from module 1 to module 3. Indeed, the second module calls the 
first module; the third module calls the second module which then calls the 
first module. A front-end program, Frontal is used to produce 3-D plots 
and contour plots of induced velocitycomputed by Velocity. All three test 
programs and Frontal use the same input file called Input. 

A. 2 First module 

The first module Location finds the critical locations in the trailing vortex 
system for a given point P and for a given rotor position V’r- Test-Location 
shows these locations in the trailing vortex system cis seen in projection in 
the tip path plane (Figure 10). Each segment is indexed by iuade^ iN and 
another index j = 1,2, respectively. There are two segments {j = 1,2) per 
spiral turn. The input data is the following and is input through the unique 
input file Input: 


• Screen: 
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0 I-PLOTTER: Indicator of plot routing (0 for screen, 1 for plot- 
ter). 

• Point position; the position of point P: 

0 rp: Radial coordinate of point P, normalized by R. 

0 Azimuth coordinate of point P in degrees. 
o zp: z coordinate of point P, normalized by R. 

• Rotor position 

0 rpfi Azimuth angle of blade of reference, in degrees. 

• Rotor characteristics: 

0 b Number of blades 

0 r„ Contracted radius, normalized by R. The trailing vortex sys- 
tem is trailed from the radial locations. 

0 N Number of spiral turns used for the trailing vortex system of 
each blade. 

• Flight conditions: 

0 /Xx Advance ratio, (F sin a/ ft R) in tip path planereference frame. 

A.3 Second Module 

The second module Wake computes the momentum theory constants for 
the given flight conditions and the ordinates of the critical locations in the 
trailing vortex system. Test-Wake shows the geometry of the distorted 
trailing vortex system as seen in perspective from downstream (Figure 11 ). 
The input data is the following and is input through the unique input file 
Input: 

• Screen: 

0 I-PLOTTER: Indicator of plot routing (0 for screen, 1 for plot- 
ter). 
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• PL0T3D: 

o PHI: Angle between i axis and vertical plane containing view 
point direction (degrees). 

o THETA: Angle between horizontal plane and view point direc- 
tion. 

0 RADIUS: distance of view point to the figure (the farther the 
view point, the most isometric perspective is obtained), (RA- 
DIUS =25 is standard)). 

• Rotor position 

0 '0,: Azimuth angle of blade of reference, in degrees. 

• Rotor characteristics: 

0 b Number of blades 

0 r„ Contracted radius, normalized by R. The trailing vortex sys- 
tem is trailed from the r„ radial locations. 

o N Number of spiral turns used for the trailing vortex system of 
each blade. 

• Flight conditions: 

0 Ct Thrust coefficient 

o fix Advance ratio, (F sina/ClR) in tip path planereference frame. 

0 Complementary ratio, {V sina/fiJR) in tip path plane refer- 
ence frame, negative in forward flight or climb. 

A. 4 Third Module 

The third module Velocity computes the induced velocity at a ^ven point 

P by the trailing vortex system and optionally by the bound vortex system. 
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A. 4.1 Test- Velocity 

Test-Velocity runs the module to be numerically tested with the debug- 
ger of the compiler. The computations are indeed traced out by extensive 
dumping throughout execution, of variable and parameter values. These 
values are then matched with hand calculations and correlation with geo- 
metrical insight available through the graphics generated by the first two 
modules. The results are given in term of A*, reduced induced inflow ratio, 

A* = Aj/A,oH 

where Aj is the induced inflow ratio: 


\i = v/ClR 


and Xifjff is the momentum value of induced inflow ratio at hover: 

XioH = — \/C'r/2 


Thus, a positive value corresponds to an induced downwash and a value 
of 1 corresponds to Xi = A^oh = —)JCt/2 = —0.056 The input data is the 
following and is input through the unique input file Input: 


• Point position; the position of point P: 


o rp: Radial coordinate of point P, normalized by R. 
o tj)p: Azimuth coordinate of point P in degrees, 
o 2 p: z coordinate of point P, normalized by R. 


• Rotor position 

; o Azimuth angle of blade of reference, in degrees. 

' • Rotor characteristics: 

o b Number of blades 

I 

I o r„ Contracted radius, normalized by R. The trailing vortex sys- 

^ tern is trailed from the r„ radial locations. 

I o cp Vortex core radius, normalized by R (cp = Ccc') 

i 

j 
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o N Number of spiral turns used for the trailing vortex system of 
each blade, 

• Flight conditions: 

o Ct Thrust coefficient 

o fix Advance ratio, (F sina/flil) in tip path planereference frame. 

o fi^ Complementary ratio, (F sina/flJ?) in tip path plane refer- 
ence frame, negative in forward flight or climb. 

• Circulation: 

o coc Zeroth harmonic coefficient for the bound circulation 

o oic First harmonic, longitudinal, for the bound circulation 

o bic First harmonic, lateral, for the bound circulation 

o d Weighting coefficient relative to remainder term contribution 
to induced downwash computation 

• Options: 

0 I AVERAGE Option control integer relative to wether or not radial 
averaging computation of induced downwash is used 

o Ar Radial length used to average the induced downwash when 
such averaging is used {Iaverage — 1)> (normalized by R). 

0 Ibound Option control integer relative to wether or not bound 
circulation induced downwash contribution is computed. 

A.4.2 Frontal 

Frontal runs the module to compute induced downwash on grid points of 
the flow field with or without averaging through the revolution of the rotor: 
A value Nmean of 1 corresponds to an induced downwash computation at 
one point of the flow field with the axis of alignment of the reference blade 
passing through that blade. It is thus the instantaneous value of the induced 
downwash seen on the reference blade. A value Nmean of 2 corresponds 
to the previous induced downwash computation averaged with another one 



A USER’S MANUAL 


44 


made after the rotor has revolved one half sector. Experimental data can 
thus be reproduced by choosing an appropriate Nmean (5 for example). 
This program produces both three dimensional (Figure 12, 13 and 14) and 
contour plots (Figure 15 etc. ) of the results. The input data is the following 
and is input through the unique input file Input: 

• File: 

0 I-INPUT: Indicator of whether an output data file is to be read 
( I-INPUT = 1 ) or new computations are to be performed ( 
I-INPUT = 0 ). 

• Screen: 

0 I-PLOTTER: Indicator of plot routing (0 for screen, 1 for plot- 
ter). 

• Contour: 

o I-CONTOUR: Indicator of plot made (0 for 3-D plot, 1 for con- 
tour plot). 

• PLOT3D: Relevant for 3-D plots only 

— ♦ 

0 PHI: Angle between i axis and vertical plane containing view 
point direction (degrees). 

0 THETA: Angle between horizontal plane and view point direc- 
tion. 

0 RADIUS: distance of view point to the figure (20 gives almost 
an isometric perspective) 

• Rotor characteristics: 

0 h Number of blades 

0 Contracted radius, normalized by R. The trailing vortex sys- 
tem is trailed from the radial locations. 

0 cr Vortex core radius, normalized by R (c/j = Ccc') 
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0 N Number of spiral turns used for the trailing vortex system of 
each blade. 

• Flight conditions: 

o Ct Thrust coefficient 

0 Advance ratio, {V sina/fli?) in tip path planereference frame. 

0 Complementary ratio, (F sina/Oi?) in tip path plane refer- 
ence frame, negative in forward flight or chmb. 

• Circulation: 

o coc Zeroth harmonic coefficient for the bound circulation 

o aic First harmonic, longitudinal, for the bound circulation 

0 bic Jb'irst harmonic, lateral, for the bound circulation 

o d Weighting coefficient relative to remainder term contribution 
to induced downwash computation 

• Options: 

o I AVERAGE Option control integer relative to wether or not radial 
averaging computation of induced downwash is used 

o Ar Ra.dial length used to average the induced downwash when 
such averaging is used {Iaverage = 1)> (normalized by R). 

o Ibound Option control integer relative to wether or not bound 
circffiation induced downwash contribution is computed. 

• Field: 

o radius-min: Minimum radius where induced downwash is com- 
puted, (normalized by R). 

o radius-max: Maximum radius where induced downwash is com- 
puted, (normalized by R). 

o azimuth-min-d: Minimum azimuth where induced downwash is 
computed, (in degrees). 
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o azimuth-max-d: Maximum azimuth where induced downwash is 
computed, (in degrees). 

• Mean: 

o N mean Number of divisions of one sector used to average the in- 
duced downwash over azimuth. (One sector is the angle between 
two successive blades, i.e. (2ir/b) 




Figure 10: First module visuailisation 
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Figure 11: Second module visualisation 
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Figure 12; Instantaneous induced velocities (3-D) 
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Figure 16: Induced downwaish field (Contour) 
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Figure 17: Radial averaging influence 
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Figure 21: Vortex core size influence 


Downstream — » 


Induced downwash field! 


o o o 


Vo# v-^ 

o o 

o \p 

-^00^1 

CV3 O 

°dl 

1! I! 

• c 

.l^D ' 

flj (d .Q 'd 1 

o 

o 

o 

o 

R: 

Q S 
o £r 

O 

o 9 


^ o 

"J 

i'.". 


So| ! 

o in o 

7 lit Wh 

■zcTa-a. 


0-® — 

'---.o.tt-'Y 

j r 

f ! aif.r. i I ^ \ ^ 


\ /^„-'«-o- 




/ \ 


\ / 
/ 




*^/M\ / I j '®’0 ^ 

jvx'A)//' y 

1l.*Y«\f . - 

li * I * \ 

I lU *11' 'v L 

,V^K \ N 

.^^4 '0»‘l ■» ^ 

\ o ?' V 

- V\ o-' \ '' ''' 

^\\v\ \ 




' 'n ''m/ 

I 


' //vH'-' 

'^'Yim 
')M V 


• \ 
> ' 


\ O/ < 

/ //'‘n 

'' N^- 

A 1 ' / H ' 

) J /j 


‘S.'" ^ J / ^ ' J 

' - / // (~ ■ 
'0*^ / // .n" 




« — apts g uioua Apv 


Figure 22: Freest ream inflow component influence 


Downstream- — » 





-dpts 8x n atreApY 


size influence (small value) 








ri «'o i‘o «o e*o o*o ro- e-o- q-o- 4*o- 6*o- ri- 

« — apis SuptraApv 


Figure 24: Remainder term size influence (large value) 
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Figure 25: Circulation zeroth harmonic component influence (small value) 
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Figure 26: Circulation zeroth harmonic component influence (large value) 
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Figure 29: Circulation first harmonic component influence on induced 

downwash field (longitudinal) 
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Figure 31: Combination of circulation first harmonic components 
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Figure 32 : Blade bound circulation influence 
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Figure 33: Final choice of circulation first harmonic components 
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Figure 34: Experimental results 
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Figure 35: Comparison computation-experiment at 0. degrees azimuth 
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Figure 36: Comparison computation-experiment at 30. degrees azimuth 
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Figure 37: Compaxison computation-experiment at 60. degrees azimuth 
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Figure 38: Comparison computation-experiment at 90. degrees azimuth 
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Figure 39: Comparison computation-experiment at 120. degrees azimuth 
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Figure 40: Comparison computation-experiment at 150. degrees azimuth 
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Figure 41: Comparison computation-experiment at 180. degrees azimuth 
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Figure 42: Comparison computation-experiment at 210. degrees azimuth 
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Figure 43: Comparison computation-experiment at 240. degrees azimuth 
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Figure 44: Comparison computation-experiment at 270. degrees cizimuth 
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Figure 45: Comparison computation-experiment at 300. degrees azimuth 
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Figure 46: Comparison computation-experiment at 330. degrees azimuth 
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Figure 47: Comparison Free- Wake and experiment at 0. degrees azimuth 
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Figure 48: Comparison Free- Wake and experiment at 30. degrees azimuth 
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Figure 49: Comparison Free- Wake and experiment at 60. degrees azimuth 
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Figure 50: Compaxison Free- Wake and experiment at 90. degrees azimuth 
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Figure 51: Comparison Free- Wake and experiment at 120. degrees aaimuth 
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Figure 52: Comparison Free- Wake and experiment, at 150. degrees azimuth 
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Figure 53: Comparison Free- Wake 
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experiment at 180. degrees azimuth 
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Figure 54: Comparison Free- Wake and experiment at 210. degrees azimuth 
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Figure 55: Comparison Free- Wake and experiment at 240. degrees azimuth 
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Figure 56: Comparison Free- Wake and 


experiment at 270. degrees azimuth 
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Figure 57: Comparison Free-WaJce and experiment at 300. degrees azimuth 
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